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Abstract 

The short parton production phase in high-energy heavy-ion coUisions is treated 
analyticahy as a nonhnear diffusion process. The initiai buiidup of the rapidity 
density distributions of produced charged hadrons within 2± 0.25 fm/c occurs 
in three sources during the colored partonic phase. In a two-step approach, the 
subsequent diffusion in pseudorapidity space during the interaction time of Tint — 
7-10 fm/c (mean duration of the colhsion) is essentially linear as expressed in the 
Relativistic Diffusion Model (RDM) which yields excellent agreement with the data 
at RHIC energies, and allows for predictions at LHC energies. Results for d-|-Au are 
discussed in detail. 
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1 Introduction 

A time-dependent analytical description of particle production from the avail- 
able relativistic energy in heavy-ion collisions is of considerable interest. In 
particular, the accurate modeling of transverse momentum and rapidity dis- 
tribution functions for produced particles is a basic requirement in attempts 
to understand the relevant partonic and hadronic physical processes. Analyt- 
ically solvable models offer transparent approaches to the problem, including 
the possibility to extrapolate to other incident energies, such as from RHIC 
energies ./s^^lQ.Q - 200 GeV to LHC, -^/s^^ 5.52 TeV. 

In this work I propose a nonequilibrium-statistical (and hence, time-dependent) 
approach that provides analytical solutions for pseudorapidity distributions of 
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produced particles. It is based on a nonlinear diffusion equation in rapidity 
space, which accounts for the explicit dependence of the diffusion coefficient 
on the rapidity density in the initial short (xp ~ 0.25 fm/c) partonic phase 
of the collision when most of the particles are produced. This is followed by 
the somewhat more extended phase of color neutralization [16], and a long- 
lasting color-neutral pre-hadronic or hadronic phase with rapid expansion of 
the system. 

Three sources for particle production are considered, two for initial rapidi- 
ties close to the beam values, and a third central source that arises mostly 
from gluon-gluon collisions. I consider the transition from the initial, highly 
nonlinear partonic phase to the subsequent, essentially linear phase. Other 
investigations such as p!7], and references to numerical simulations therein, 
corroborate the short duration of the parton-production phase, and the long 
duration of the subsequent recombination (pre-hadronic) and hadronic phase. 
According to parton-cascade models such as p8], the pre-hadrons are color 
singlets that are generated from quark and gluon recombination in a statisti- 
cal coalescence process. They decay into the final hadrons according to their 
relative phase-space weights. Many numerical approaches to the problem use 
string models in the initial phase, and final-state hadronic collisions, but often 
hadronic rescattering is not considered. 

The second phase lasts about Tint — 7-10 fm/c (depending on the system, the 
incident energy, and the centrality). Phenomenologically, it is well accounted 
for in a Relativistic Diffusion Model (RDM) based on a linear Fokker-Planck 
equation with constant diffusion coefficient, plus fast collective expansion. 
The RDM had been developed some time ago p!|2|l3lH|5] and compared in 
detail with data on net proton rapidity distributions at SIS, AGS, SPS, and 
RHIC energies, and with produced-particle distributions. For pseudorapidity 
distributions of produced charged hadrons, a x^-optimization of the Jacobi- 
transformed analytical solutions yields very precise agreement with the avail- 
able RHIC data provided the midrapidity source for particle production is 
taken into account. 

For net-proton rapidity distributions at RHIC energies, there have been several 
investigations of nonlinear effects within the diffusion approach. A nonlinearity 
in the drift coefficient that secures the correct Maxwell-Boltzmann equilibrium 
limit for t —>■ oo has been investigated in [8], but the deviations from linearity 
are small. 

The strong effect of a nonlinearity in the diffusion coefficient that also persists 
over the full interaction time of typically 10 fm/c [Tip5] as a consequence of 
the introduction of non-extensive statistics [10] had been investigated in [8|39] 
for net protons in heavy-ion collisions, and in [U] for produced particles in pp- 
coUisions. It seems that this is a way to account for the collective expansion 
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of the system without considering it exphcitly. If one includes an exphcit 
treatment of collective expansion |3D], however, there is no need to introduce 
non-extensive statistics when comparing the RDM-results with data, the linear 
evolution after the parton-production phase yields excellent results for both 
net baryons, and produced particles. 

The linear diffusion model had also been proposed in order to calculate and 
predict transverse energy distributions of hadrons |6], and more recently to 
calculate transverse momentum distributions of identified hadrons (neutral 
pions, as well as negative pions and kaons) in a nonequilibrium-statistical 
approach including radial flow [7j. 

Equilibrium-statistical models [TT]|12|I13] account in remarkable detail for rel- 
ative production rates of produced particles at central rapidity with only the 
temperature and the chemical potential as parameters. Due to the lack of 
time dependence - and consequently, of nonequilibrium-statistical effects - a 
relevant ingredient is, however, missing if one aims at the precise modeling of 
distribution functions. 

The gradual approach of the system towards statistical equilibrium in the 
course of a relativistic heavy-ion collisions is presented in a schematic analyt- 
ical model in this work. I start with the nonlinear diffusion model as expected 
to be valid during the parton-production phase in Section |2} followed by a 
consideration of the so-called source solution of the nonlinear problem in Sec- 
tion |3| For sufficiently large times, the initial power-law behaviour during 
the nonlinear parton-production phase is superseded by the essentially linear 
diffusion (RDM-) phase in rapidity space which produces gaussian tails in 
pseudorapidity space. Section |4j 

The late-stage time evolution is discussed in a comparison with RHIC data 
for the asymmetric d -|- Au- system, which indeed show the gaussian tails. In 
an asymmetric system like d-l-Au, the nonequilibrium effects are visible more 
directly than in case of symmetric systems such as Au -|- Au. The interaction 
ceases long before statistical equihbrium with respect to the variable rapidity 
is reached. The conclusions are drawn in Section [5l 



2 The nonlinear diffusion equation 

The origin of diffusion during and after particle production in a heavy-ion 
reaction at relativistic energies is found in momentum space, through random 
momentum kicks of the produced particles - partons in the soft-gluon field in 
the early stage, prehadrons and hadrons in later stages. Diffusion in coordinate 
space appears as a secondary effect. 



3 



The corresponding fluctuations can be seen in rapidity and pseudorapidity 
distributions of produced particles [TII2]|S|^|f5] . as well as in transverse energy 
and momentum distributions pfT] . To provide an analytical treatment of the 
problem in both early and late stages, I confine the present work to rapidity 
space, with the lorentz-invariant rapidity y = 0.5 ■ ln{{E + p)/{E — p)), and a 
subsequent Jacobian transformation to pseudorapidity r] that is required for 
the comparison to the available data for produced charged hadrons. 

To incorporate the early parton-production phase into the relativistic diffusion 
model [USE], a dependence of the diffusion coefficient on the initially very 
high rapidity density R{y, t) has to be considered, Dy Dy{R), such that the 
linear transport equation in rapidity space that I investigated in [1] is replaced 
by 

^^R{y,t) = -Vy[j{y)R{y,t)]+VyDy{R{y,t))VyR{y,t). (1) 



The drift term J{y) governs the gradual approach of the mean values towards 
statistical equilibrium. The diffusion coefficient Dy{R) depends on the rapidity 
density and hence, the equation is generally highly nonlinear. It is therefore 
expected to account not only for the long-lasting, essentially linear diffusive 
phase as in p!|2|3] . but also for the partonic initial phase of high rapidity 
density. In this short-lived phase the major part of particle production with 
rapidly rising norm of the distribution function takes place. The rising norm 
is phenomenologically accounted for in this work by letting the integration 
constant in Eq.Q depend on particle number. 

In case of the linear RDM [li6 \ with Dy[R) = Dy=const., I had assumed an 
instant production of the particles in the three sources, and subsequent diffu- 
sion in ?/— space during the interaction time. This initial condition is exactly 
fulfilled only for net baryons [T]. However, for produced charged hadrons, it 
also yields extremely precise results when compared in detail to the 

data. An explicit treatment of the short nonlinear parton production phase 
with a strong dependence of the diffusion coefficient on the distribution func- 
tion should therefore preserve the model features of the subsequent, essentially 
linear diffusive phase. 

To account for the strong correlation between diffusion coefficient and rapidity 
density distribution in the initial high-density particle production phase, I 
propose a dependence on a power k, of the rapidity density according to 

Dy[R{y,t)]=Dl.Riy,tr (2) 



with the rapidity diffusion constant in the particle production phase. For 
certain critical exponents k,, analytical solutions of the diffusive part of the 
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transport equation can be obtained. 



In a two-step approach, the subsequent - probably mostly pre-hadronic and 
hadronic - evolution in pseudorapidity space during the interaction time of 
Tint — 7-10 fm/c (mean duration of the collision; see [H] ) is accounted for 
in the Relativistic Diffusion Model [1] (RDM, k = 0) with a linear drift term 
Jiy) = iVeq — y)l^y govemed by the rapidity relaxation time Ty and the equi- 
hbrium value of the rapidity Deq- The diffusion term is here oc Dy^R{y,t). 

The diffusion constant Dy in this phase is significantly smaller than the value 
of in the short production phase with a large number of (partonic) degrees 
of freedom. The linear model with instant particle production yields excellent 
agreement with d+Au, Cu+Cu and Au+Au data at RHIC energies, including 
the detailed centrality dependence [3]. 

The initial short, highly nonlinear phase of parton production occurs within 
Tp ~ 0.25 fm/c [1^ in two beam-like sources, and a central source in rapidity 
space. Since the time scale for particle production in all three sources is faster 
than the one for the nonlinear diffusion, it turns out that the particle content 
in the power-law tails remains small: there is little spread of the distribution 
function in rapidity space in this initial phase. 

The mathematical treatment of the initial nonlinear phase is confined to the 
central source as an example. The beam-like sources are then dealt with in an 
analogous way, but with the freedom to choose different diffusion coefficients 
because the production mechanisms in the valence-quark dominated beam- 
like regions of rapidity space are different from the central region with few 
valence quarks at RHIC or LHC energies. Kinematic constraints in the beam- 
like sources expected at high absolute values of rapidity are not considered 
here. 

During the short production phase, the drift in ?/— space is not yet pronounced, 
and I therefore treat here only the diffusive part of the nonlinear transport 
equation, R{y, t) Piv, t) with 



The solution of this equation at the end of the nonlinear production phase 
{t = Tp) can then be used as initial condition for the subsequent linear diffusive 
time evolution treated in the Relativistic Diffusion Model (RDM) [Tll2ll3llill5] . 



dt 



P{y,t) = DP^VyP{y,trVyP{y,t). 



(3) 



Ro{y,0)=Piy,Tp). 



(4) 
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With t* = t ■ the nonhnear diffusion equation becomes 

^P{y.t*) = VyP{y,erVyP{y,f). (5) 

This equation has been extensively studied in many diverse areas of science, 
and a large amount of mathematical literature exists, P^20f21f22ll23f24f25ll26|l2' 
and references therein. It has mostly been considered for positive values of k 
such as K = 1 for thin saturated regions in porous media, k > 1 for the per- 
colation of gas through porous media, k = 3 for thin films spreading under 
gravity, and k = 6 for radiative heat transfer by Marshak waves. Many prob- 
lems are dealt with in only one (spatial) dimension, analogous to the present 
work which is confined to one (momentum-like) dimension. It has also been 
established that the number of exact solutions is limited. Solutions for sev- 
eral power-law diffusivities with negative k are known [25j, in particular, for 
K = -1/2, -1, -4/3, -3/2 and -2. 

The physically most interesting behaviour occurs in the region of small dif- 
fusion coefficients D(R), where a moving boundary may exist. The behaviour 
of solutions for positive and negative values of k is distinctively different. De- 
pending on the specific values for the constants of integration (see below), 
a free boundary may occur for n = 1 that has a finite gradient and moves 
with finite velocity, similarly for other positive values of k, but for k > 1, the 
gradient at the boundary becomes infinite. 

For K = —1 and some other negative values, however, there is an instanta- 
neous spread without a free boundary, as was treated by Pattle [12], Pert 
[22] Hill [25] and others for instantaneous heat deposition in a medium with 
concentration-dependent diffusion coefficient. This situation is in some respect 
analogous to the initial parton production from the available energy in a rela- 
tivistic heavy-ion collision. Here the initial rapidity density distribution of the 
created partons should not have a free boundary in rapidity space, except for 
kinematical constraints. 

Hence I investigate solutions for negative values of k emphasizing k = —1. An 
exact solution is not only useful to test the accuracy of numerical results, but 
it is also important to understand and describe the physical behaviour of the 
system. In particular, the analytical solution at the end of the initial nonlinear 
phase may then be used as initial condition for the subsequent, essentially 
linear diffusion process in y— space which can be modeled analytically within 
the given simple but successful RDM-framework. 

The majority of known exact solutions of the nonlinear problem are so-called 
similarity solutions [19f2UII21f22f23ll24f25f26j : With an assumed functional form 
of the solution, the partial differential equation reduces to an ordinary differ- 
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ential equation, or to a partial differential equation of lower order, which can 
then be integrated in closed form under certain conditions. The similarity 
solution of ([5]) for t* — t is written as 

P{y,t)=y''/^^'+'^m (6) 



with 

e = (7) 

and A is an arbitrary constant with A 7^ —1. Inserting this ansatz into the 
nonlinear differential equation ([s]) yields first integrals for two values of A (for 
K = —1 and —2, only one value of A) 

Xi = -k/{k + 2) (8) 

X2 = -k/{k + 1) (9) 

which have been given by Hill and Hill (1990) [23] as 
A'^*' 2$"+^ 2$ 

+ . , OA. = Cl (10) 



{2k + 3)$"+^ $ 
^ («:+ 1)2^2 +(2«:+l)2- ^ ^ ^ 



with $' = d^/d^. For vanishing constants of integration Ci = 0, the so- 
called "source" and "dipole" solutions arise from these first integrals. For 
finite Ci 7^ 0, a number of exact solutions for special values of k, have also 
been derived [25]. Here I investigate the source solution for Ci = 

$ = (C2^W(«+2) _ (12) 



This solution is not defined for k = —2, and becomes singular for n < —2. 
Hence, a solution without a free boundary as required for the initial rapidity 
diffusion problem can in principle only occur for —2 < n < 0. The desired 
solution for K = — 1 (and hence, A = Ai = 1) obeys the first-order partial 
differential equation 

$'2 

^-^ + 2$ = (13) 
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and the solution of the nonhnear diffusion equation 



(14) 



with ^ = {y/ty^"^ becomes 

P{y,t) = [C2yr' + yem-' (15) 

where C2 is the constant of integration. 



3 The source solution in particle production 

The solution of the nonlinear diffusion problem in a high-density phase such 
as during parton production in rapidity space has thus been reduced to 

Piy,t) = [C,t + yy{2t)]-\ (16) 

To account for the increasing norm of the rapidity distribution function during 
the rapid parton production process, I let the integration constant C2 depend 
on the particle number. Since the particle number is likely to increase expo- 
nentially with time during the first collision phase where t = t < ~ 0.25 
fm/c, I take the particle-number dependence into account phenomenologically 
by choosing the integration constant in the denominator as 

C2=exp(-t). (17) 



In the present two-step model, all particles are assumed to be produced until 
t = Tp where the nonlinear production phase turns into an (essentially linear) 
diffusion process with k = 0. That is, for C2 = exp(— Xp) the norm of the 
solution reaches its maximum value 

J P{y,t = Tp)dy = l. (18) 



The corresponding solutions ( 16 ) of the nonlinear diffusion problem are shown 
in Fig. 1 for the central source of parton production in d+Au at y/sNN= 
200 GeV. Results of the rapidity distribution functions for produced particles 
multiplied with the number of charged hadrons produced in the central source 
{N^f^ = 22, 13]) are shown for four values oit* = t- Dp=U, 14.5, 15, 16.4, with 
the norm of P(y,t) reaching / P{y,t)dy=l for t=rp=0.25 fm/c at t*=16.4. 
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Fig. 1. Time evolution of rapidity distributions for produced charged particles from 
minimum-bias d + Au collisions at -^/snn = 200 GeV in the initial nonlinear phase 
for the midrapidity source. The norm of the distribution is rising with increasing 
t* = Dy-t until at t = Tp all the particles have been produced and the linear diffusion 
process in y— space starts. The dashed curve is a Gaussian with the same particle 
content N^f^ = 22 as the power-law solution of the nonlinear problem at t = Tp. 



The distribution functions are seen to rise strongly with increasing t* in a 
narrow midrapidity region due to the fast increase of the particle number, 
with an only moderate increase in the power-law tails where the particles 
with higher rapidities are created: these tails arc already present at very short 
times. Physically, the partons with the highest rapidity values are created 
already at the shortest times. 



The rapidity diffusion coefficient in the particle production phase is thus 
= 16.4/rp ~ (16.4/0.25)c/fm ~ 66c/fm, which is significantly larger than 
the diffusion coefficient Dy in the subsequent long-lasting (~ 7 fm/c) linear 
diffusion phase. This reflects the larger number of degrees of freedom in the 
initial phase which is mainly partonic, and the higher density of particles. 



The dashed curve in Fig.l is a Gaussian that arises from a linear time evolution 
{k=0) and 5— function initial conditions with instant particle production at 
t = 0. It has the same particle content (22 charged hadrons) as the nonlinear 
solution at t* = TpDP = 16.4. 
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For the beam-like sources at initial rapidities yi^2 = TUmax, the solutions of 
the nonlinear problem for k = — 1 are accordingly {t* = t ■ -D^i 2 ~^ t) 

Pi,2{y,t) = [C2t + {y±yma.?im]-\ (19) 



Here the diffusion coefficients in the particle production phase D^i^2 likely 
to differ from in the midrapidity region because the microscopic processes 
during particle production are substantially different, with the diffusion co- 
efficient at midrapidity mostly due to gluon-gluon collisions. The analytical 
solutions at the end of the initial phase t = in the three sources are dis- 
played in the upper part of Fig. 2 for minimum-bias d+Au at ^JsNN= 200 
GeV. 

The corresponding number of charged hadrons [3] created in the Au-like source 
is Nly^=bb, in the d-like source iV2^=14, and in the central source N^f^=22, with 
a total of 91 produced charged hadrons in minimum-bias collisions. These 
particle numbers have been determined from a detailed comparison of the 
subsequent linear diffusion phase with data, see Section |4} In the lower part 
of Fig. 2, Gaussians with the same particle-number content are displayed, as 
they arise from a linear diffusive time evolution (k = 0) with 5— function initial 
conditions, see next section. 

At t = Tp the power-law solution of the initial nonlinear phase in the midra- 
pidity source can be expressed as 

P3{y,t) = ^^^ (20) 



with 



C = 2TpDl 



(21) 
(22) 



for the central source, and analogously -Pi,2(y, t) for the beam-like sources with 
y ^y± ymax and Dp D^^^- 



As was shown in [Tp|3Kll5] . the linear diffusive evolution after the initial short 
parton production phase is in very good agreement with the available data. 



Hence the power-law result (20 ) of the first phase can be used as an initial con- 



dition for the second, essentially linear phase that is reconsidered as described 



by Eq. (31) in the next section. With this initial condition, the solution of the 
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Fig. 2. Rapidity distributions for produced charged particles from minimum-bias d 
+ Au collisions at y/s^N = 200 GeV at the end of the initial nonlinear phase (top) 
for the three sources, with particle contents N^^ = 55,N^f^ = 14, and N^^^ = 22. 
In the bottom frame, corresponding Gaussians are displayed (see also Figjs]). They 
arise from a linear evolution with (5— function initial conditions. 
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linear diffusion problem for the central source becomes |28] 



R{y,t) 



C 
2a 



\ 



2tt 



0", 



3? 



exp 



{y + ivy 



erfc( 



ly — V ^ 



(23) 



with the variance 

2 



a^(t)=Z}^r,[l-exp(-2t/r,)] 



(24) 



and 



V = aexp(— t/r. 



(25) 



For t 



3? 



^> oo we have 
erfc( ^ 



(26) 



such that the gaussian limit is attained 



R{y,t 



oo 



\ 



27ra2 



exp 



(27) 



which is a solution of a linear Fokker-Planck equation (31) for t ^ oo with 



(5— function initial conditions. For gaussian initial conditions with an initial 
variance cTq, one obtains a gaussian solution with o"^(t) —>■ ay{t)-\-o'lexp{—2t/Ty) 



and hence, Eq.(27) results for t — > oo as well. 



To connect the diffusion approach with data, the linear Relativistic Diffusion 
Model (RDM) for the second, long-lasting diffusion phase in pseudorapidity 
space [Tpll3llill5] is reviewed in the next section. 



4 Linear diffusion phase 



Since the initial power-law behaviour is superseded by the Gaussian evolution 
at sufficiently large times, the evolution is started here at t=0 with 5— function 
or gaussian initial conditions to illustrate the outcome of the three-sources 
model for large times and in particular, to compare to data. 
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The situation at moderate times with 5— function initial conditions is shown in 
the lower part of Fig. 2, with separate Gaussians in the three sources in rapidity 
space which have the same particle-number content as the initial power-law 
solutions of the nonlinear particle-production problem. The subsequent time 
evolution of these three sources leads to agreement with the experimental data 
at the interaction time t = Tint, and to statistical equilibrium for t ^ oo. 

The nonequilibrium-statistical description of this evolution is based on an 
essentially linear diffusion equation which is briefly reviewed in this section. 
We have used a Fokker-Planck equation (Uhlenbeck-Ornstein [22] version with 
K = 0) [Tll2|3pi)5] for the distribution function R{y, t) for produced charged 
hadrons in rapidity space 

^^R{y,t) = -Vy[j{y)R{y,t)] + DyVlR{y,t). (28) 



The drift is now taken into account since we look at the large-time behaviour, 
and the drift function J{y) determines the speed of the statistical equilibration 
in y-space. In order to attain the Boltzmann distribution for large times, the 
drift term must have the form pi39] 

J{y) oc m_L sinh(?/) oc p\\ (29) 



with the transverse mass m_|_ = y {m? +p±), and the longitudinal momentum 
This introduces another nonlinearity into the problem, which prohibits an 
analytical solution of the diffusion equation. Such an analytical solution pj is, 
however, possible by linearising the drift function in a relaxation ansatz 

J{y) = {ye,-y)/Ty (30) 



with the rapidity relaxation time Ty, and the equilibrium value of the rapidity 
yeq that is calculated from energy and momentum conservation in the system of 
participants. The deviations of the solutions for nonlinear and linear versions 
of the drift are not pronounced [32] and hence, I have used the analytical 
solutions of the linear problem for the components Rk{y,t) (fc=l,2,3) of the 
distribution function 



d_ 
di 



Rk{y,t) 



-Tydy 



{y-Veq) ■ Rkiy,t) 



D';^,Ri^{y,t). 



(31) 



The diagonal components Dy of the rapidity diffusion tensor contain the mi- 
croscopic physics in the respective beam-like {k = 1,2) and central {k = 3) 
regions. They account for the statistical broadening of the distribution func- 
tions. To connect with data, one has to consider the additional broadening 
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due to longitudinal collective expansion that leads to a larger (effective) value 
of Dy [30] than what is calculated [T] from the dissipation-fluctuation theorem 
(Einstein relation). 

As discussed above, the initial conditions in the linear phase are taken as 
Ri,2{y,t = 0) = 6{y ± Umax) with the maximum rapidity i/max = 5.36 at the 
highest RHIC energy of y/s^N = 200 GeV (beam rapidities are yi^2 = TUmax), 
and R-iiy.t = 0) = Sijj). A midrapidity gluon-dominated symmetric source 
had also been proposed by Bialas and Czyz [30] . 

This initial condition for the midrapidity source in the linear case corresponds 
to initial particle production without any longitudinal motion, independently 
of the mass of the collision partners: the third source is created at |/ = 0, 
and the drift towards the equihbrium value y = yeq, as well as the rapid 
collective expansion, sets in subsequently. (In contrast, the nonlinear model 
as discussed in the previous section produces power-law tails at short times, 
which are superseded by the gaussian tails of the linear evolution only at later 
times). 

The mean values in the three sources have the time dependence 

< 1/1,2 (i) >= 1/eJl - exp(-t/ry)] =F y-max exp {-t/Ty) (32) 



for the sources (1) and (2), and 

< y^it) >= l/ejl - exp(-t/r,)] (33) 



for the moving central source. The three mean values reach the equilibrium 
limit for time to inflnity. In our previous RDM-calculation [3j with slightly 
different initial condition, the mean value of the central source was at the 
equihbrium limit < ysit) >= y^q independently of time, thus assuming instant 
equilibration in this source regarding the mean values. The variances erf 2 sit) 



are as in Eq.(24), with Dy ^ Dp. 



It turns out that for d+Au at the highest RHIC energy, one can not de- 
termine from a comparison with the data which of the two possibilities for 
the initial conditions of the central source is more realistic because the is 
nearly identical in both cases. The subsequent diffusion-model time evolution 
in pseudorapidity space is followed up to the interaction time Tint, when the 
produced charged hadrons cease to interact strongly. 

The quotient Tint/Ty is determined from the minimum with respect to 
the data, simultaneously with the minimization of the other free parameters 
- namely, the variances of the three partial distribution functions, and the 
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Table 1 

Produced charged hadrons in minimum-bias d + Au collisions at ^snn = 200 
GeV, yi^2 = T 5.36 in the linear Relativistic Diffusion Model. The equilibrium 
value of the rapidity in the RDM is yeq, the time parameter (see text) is p, the 
corresponding value of interaction time over relaxation time is Tint/Ty, the variance 
of the central source in y— space is cTo. The number of produced charged particles 

12 

IS NX for the sources 1 and 2 and N^^^ for the central source, the percentage of 
charged particles produced in the midrapidity source is nj?^, and x^/d-o.f. is the 
result of the minimization [3] per number of degrees of freedom. 
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number of particles produced in the central source. In this nonequilibrium- 
statistical approach, the equilibrium value of the rapidity and its dependence 
on centrality is calculated from energy and momentum conservation in the 
system of participants as 



(34) 



with the transverse masses < rnj2{b) >= y("^i,2(^)+ < Pt >'^), and masses 
nii,2(^) of the "target" (Au)- and "projectile" (d)-participants that depend on 
the impact parameter b. The average numbers of participants < Ni^2{b) > 
from the Glauber calculations reported in [31] for minimum bias d + Au at 
the highest RHIC energy are < A^^i >=6.6, < N2 >=1.7, which we had also 
used in [3]. 

The average numbers of charged particles in the target- and projectile-like 

1 2 

regions N^f^ are proportional to the respective numbers of participants Ni^2, 

- ^^'^ iN, + N2) ^^^^ 



with the constraint A^*^* = A'^ + ^ch + ^ch- Here the total number of charged 
particles A^*^* is determined from the data. The average number of charged 
particles in the equilibrium source A"^^ is a free parameter that is optimized 
together with the variances and Tint/Ty in a x^"fit of the data using the CERN 
minuit-code. The results are summarized in table 1. 

The FPE is solved analytically as outlined in [1], and the solutions are con- 
verted to pseudorapidity space. This conversion is required because particle 
identification is not available. The relation between scattering angle 6 and 
pseudorapidity r] is i] = — ln[tan(6'/2)]. Here 9 is measured relative to the di- 
rection of the deuteron beam. Hence, particles that move in the direction of 
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the gold beam have negative, particles that move in the deuteron direction 
have positive pseudorapidities. The conversion from y— to rj— space of the 
rapidity density 

dN pdN ..dN 



is performed through the Jacobian 

j(r?, (m)/(pr)) = cosh(r7) ■ [1 + {{m) / (j^t)? + smh\r,)]-^'\ (37) 



The average mass < m > of produced charged hadrons in the central region 
is approximated by the pion mass m^r, and a mean transverse momentum 
<Pt > = 0.4 GeV/c is used Due to the conversion, the partial distribution 
functions are different from Gaussians. The charged-particle distribution in 
rapidity space is obtained as incoherent superposition of nonequilibrium and 



local equilibrium solutions of (31) 

dNch{y,t = Tint) 



NlhRi {y, Tint) + Nl^R2 {y, nnt) + N^hRs (y, nnt) (38) 



with the interaction time Tint (total integration time of the differential equa- 
tion). The integration is stopped at the value of Tint/Ty that produces the 
minimum with respect to the data and hence, the explicit value of Tint is 
not needed as an input. The resulting values for Tint/Ty are given in table [l] 
together with the widths of the central distributions, and the particle numbers 
in the three sources. 

The time evolution is shown together with the comparison to PHOBOS minimum- 
bias data [21] in Fig|3] It is evident that the two beam-like distribution func- 
tions move towards smaller absolute pseudorapidities as time increases, and 
reach agreement with the data at p=0.54. Here the time evolution parameter 
p in the numerical calculation is defined as p] 

p= 1 -exp(-t/rj,). (39) 



The minimum-bias result also shows the asymmetric shape of the distribution 
function, which is very well reproduced in the diffusion calculation. At larger 
values of the time evolution parameter p, all three subdistributions tend to 
become symmetric in y with respect to the equilibrium value yeq, indicating 
the approach to thermal equilibrium. At p=0.999, the equilibrium state is 

^ There is a difference of a factor of two in the exponent as compared to the 
definition of p used in [3], which causes different t/Ty values for given p. 
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Fig. 3. Time evolution of pseudorapidity distributions for produced charged particles 
from minimum-bias d + Au collisions at ^/sNN = 200 GeV in the linear diffusion 
model as outlined in |3|4] . The d-like source is shaded to illustrate the movement 
in r/— space towards equilibrium. Dash-dotted curves show the slightly moving, glu- 
on-dominated midrapidity source for hadron production. Results for five time steps 
(p-values, cf. text) are displayed. Agreement with the data [H] is reached at p=0.54. 
Statistical equilibrium centered at r/gg would be achieved at later times. Reviewed 
from |4j. 
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already closely approached. The slight asymmetry is due to the conversion 
from rapidity- to pseudorapidity space which tends to produce a dip at 77 = 0. 
For time to infinity, statistical equilibrium in pseudorapidity space would be 
reached. 

We have shown in [3] that with this linear RDM approach, the centrality 
dependence of the measured pseudorapidity distributions [31] from central to 
very peripheral collisions can also be modeled in considerable detail, Fig|4j For 
peripheral collisions, the asymmetry of the overall distribution is not yet pro- 
nounced because here the d- and the Au-like partial distributions are similar 
in size due to the small number of participants. 

Towards more central collisions, the number of gold participants rises, and 
the corresponding partial distribution of produced particles becomes more 
important. In addition, the distributions drift towards the equilibrium value. 
Both effects produce the asymmetric shape, which is also seen in minimum- 
bias. The tails of the distribution functions are of gaussian shape in perfect 
agreement with the data. This shows that the power-law tails of the initial 
phase have not survived the time evolution, as is confirmed when the result 
of the initial phase is used explicitly as an initial condition for the linear 
evolution. 

It is interesting to compare the behaviour of the rapidity or pseudorapidity 
distribution functions with results from different approaches to the problem 
such as saturation models [32)133035] . viscous hydrodynamics [37], or ideal hy- 
drodynamics [38] . 

Calculations within the framework of the Parton Saturation Model not only 
predict the midrapidity value, but also the full rapidity distribution function 
(at RHIC energies, and also at LHC) [3SE1]- These calculations are based on 
a classical effective theory that describes the gluon distribution in large nuclei 
at high energies where saturation might occur at a critical momentum scale, 
to form a Color Glass Condensate (CGC) [35] . 

This assumption has a clear and reasonable physical basis and yields good re- 
sults for pseudorapidity distribution functions of produced charged hadrons at 
the available RHIC energies. Problems may be expected for net-proton rapid- 
ity distributions since protons and antiprotons are produced in equal amounts 
from the CGC. At LHC energies, the overall pseudorapidity distribution from 
the CGC as obtained with the assumption of a constant for strong coupling 
is slightly narrower than the corresponding diffusion- model prediction [5]. 

Additional consideration of a running coupling gives a midrapidity value that 
is of the order of 10% smaller; another uncertainty arises from the extrapola- 
tion of the saturation scale to LHC energies. Various predictions for central 
rapidity densities and pseudorapidity distributions at RHIC and LHC ener- 
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r 




Fig. 4. Calculated pseudorapidity distributions of produced charged particles in d 
+ Au collisions at ^/snn = 200 GeV for five different centralities, as outlined in 
|3|4] ■ Central collisions are shown in the top frame, peripheral at the bottom. The 
linear diffusion-model (RDM) results for three sources (d-like source shaded) in 
r/— space are shown together with their incoherent sums as minimizations at 
each centrality cut (c.c). The time variable is p (see text). The initial conditions for 
the central source are slightly different from [3], see text. Data are from PHOBOS 
|31j . Reviewed from [1]. 
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gies had been summarized e.g. in [36], where also the differences among the 
existing models - including hydrodynamical and pQCD approaches and their 
numerical implementations - had been discussed. 

In relation to these approaches, the analytical diffusion model provides good 
results when compared in detail to the experimental distribution functions at 
RHIC energies, in particular, in the tails of the distributions. To provide a 
microscopic foundation, however, a derivation of the diffusion coefficients in 
the three sources would be required. Due to the valence-quark dominance in 
the beam-like sources as opposed to the mainly gluonic midrapidity source, 
the diffusion coefficients may turn out to be substantially different in the three 
sources. 



5 Conclusion 

I have presented a nonlinear Relativistic Diffusion Model that includes an 
explicit analytical treatment of the initial parton-production phase in rapidity 
space. As a consequence of the high rapidity density at short times t ~ 0.25 
fm/c, the rapidity diffusion coefficient depends on the distribution function, 
such that the problem is highly nonlinear in the initial phase. 

For a power-law dependence of the diffusion coefficient on the distribution 
function with an exponent k, the mathematical technique of similarity solu- 
tions that has been refined in many previous works proves to be useful in the 
present physical context. In particular, I have investigated the so-called source 
solution in rapidity space. 

An adequate solution to the nonlinear problem should not have a free bound- 
ary in rapidity space that moves with finite velocity as is the case for positive 
values of k, but in parton production there should be an instantaneous spread 
in I/— space without a free boundary. This corresponds to the case k = — 1, 
which can be solved analytically using the technique of similarity solutions. 

During particle production in three sources, the norm of the distribution func- 
tion increases, which I have considered phenomenologically by letting the inte- 
gration constant depend on particle number. Since particle production is very 
rapid - exponentially in time -, this increase of the norm of the distribution 
function turns out to be faster than the spread of the distribution function in 
rapidity space due to nonlinear diffusion and hence, the power-law tails of the 
distribution function remain small during the parton-production phase. 

The result of the parton-production phase is then used as initial condition for 
the later (pre-hadronic and hadronic) stage of the collision, which is treated 
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here in the hnear relativistic diffusion model. The hnear evolution washes out 
the initial power-law tails of the distribution function, which become gaussian. 



With the proper Jacobi transformation to pseudorapidity space, this approach 
yields very precise agreement with charged-hadron data for both asymmetric 



systems (d+Au), and symmetric systems such as Cu+Cu and Au+Au [2ll5] . 
It is also particularly suitable for predictions at LHC energies of 5.5 TeV for 
Pb+Pb. 

The second collision phase lasts about 7-10 fm/c depending on the system, 
the incident energy, and the centrality. Due to the schematic treatment that 
is based on a linear partial differential equation, particle number is conserved 
in this phase, which appears as a reasonable physical assumption even though 
it is certainly not strictly valid. 

For time to infinity, the evolution of the distribution function proceeds to 
statistical equilibrium with respect to the variable rapidity or pseudorapidity. 
Comparing the data with this time evolution as modeled by the solutions of 
the linear problem, it is evident that at the time of the measurement - when 
strong interaction ceases - the system is still far from statistical equilibrium. 
This underlines the view that relativistic heavy-ion collisions are very suitable 
to observe strongly interacting many-body systems with a large amount of 
particle production on their way to statistical equilbrium. 

In this work I have not considered the connection between the diffusion ap- 
proach and QCD. It is obvious that the forward and backward sources for pro- 
duced particles are related to the valence quarks, whereas the central source 
is essentially due to gluon-gluon collisions. An actual microscopic calculation 
of the three sources emphasizing their relative size (number of produced par- 
ticles) is therefore of interest ^41] . 
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